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The differential formulation of the Landau-Fokker-Planck collision integral is developed for 
the case of relativistic electromagnetic interactions. 
PACS numbers 52.25.Dg, 52.20.-j, 52.60.+h 

Kinetic theory is founded upon the Boltzmann equation, which is a conservation equa- 
tion for the phase-space distribution function of each species in an ensemble of interacting 
particles. For the case of Coulomb interactions, Landau^ expressed the collision term in the 
Fokker-Planck form. This mixed integro-differential representation was extended to rela- 
tivistic electromagnetic interactions by Beliaev and Budker.^ For the nonrelativistic case, 
it was shown by Rosenbluth et al.'^ and by Trubnikov* that the integrals appearing in the 
collision term can be expressed in terms of the solution of a pair of differential equations. 
The present work extends that formulation to the relativistic collision integral. Using an 
expansion in spherical harmonics the relativistic differential formulation is then applied to 
calculate the scattering and slowing down of fast particles in a relativistic equilibrium back- 
ground plasma. Our work is relevant to the study of high temperature plasma in fusion 
energy research and in astrophysics. 

In the work of Landau^ and that of Beliaev and Budker,^'^ the collision term that 
occurs on the right-hand side of the Boltzmann equation for species a and describes the 
effect of collisions with species b is written in the Fokker-Planck form. 

Cab — -TT- ■ {Dab ' 'Pabfa), (1) 

on ou 

in which the coefficients Dab and Yab are defined by 

Dab{n) ^^^\og Kab I (y(u,u')/b(u')d^u', 
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Here, fa and ft are the distribution functions for the two species, u is the ratio of momentum 
to species mass, qa and qt are the species charge, nia and mb are the species mass, eo is 
the vacuum dielectric permittivity, and logAaf, is the Coulomb logarithm. The kernel U 
is specified below. This form of the collision operator is only approximate because of the 
introduction of cutoffs in the collision integral. More accurate operators that take into 
account Debye shielding at large impact parameters and large-angle scattering and quantum 
effects at small impact parameters have been derived. ^'^ The purpose of this letter is to 
present a differential formulation for the integral transforms that occur in Eqs. (2). To avoid 
unnecessary clutter we discard the factor that depends only on the species properties, drop 
the species subscript, and consider the transforms 

D(u) = ^y'y(u,u')/(u')dV, (3a) 
F(u) ^-l-JI^A.. t;(u,u'))/(u')d^u'. (3b) 
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For guidance, let us recall briefly the nonrelativistic theory.^'^ In that case the 
moment um-to-mass ratios u and u' reduce to the velocities v and v', and the collision 
kernel is the one given by Landau/ U = (|sp/ — ss)/|s|^, where s = v — v'. It may be seen 
that U = 9^|s|/9v9v and {d/dv')-U = —2d\s\~-^/dv. To obtain the differential formulation, 
these representations arc inserted into Eqs. (3), and the differentiation with respect to v is 
moved outside the integration over v'. Defining the potentials h{v) = — (I/Stt) / \s\f (Pv' 
and g{v) = -(l/47r) J\s\-'^f(fV, we have D = -d'^h/dwdw and F = -dg/dw. Further- 
more, from A|s| — 2|s|~^ and A|s|~^ = —Ati8{s) it follows that h and g obey the equations 
A/i = g and Agi = /. (A denotes the Laplacian with respect to the variable v.) These 
equations provide the differential formulation of the collision term in the nonrelativistic 
case. 

The Landau collision kernel was obtained in a semi-relativistic fashion, assuming 
Coulomb collisions and relativistic particle kinematics. It is a good approximation to the 
fully relativistic kernel given below provided that |v.v'| <C c? , which is true when one of the 
colliding particles is nonrelativistic. However, the reduction of the collision integral to the 
differential form of Rosenbluth and Trubnikov relies on the stronger assumptions |v|^ <C 
and |v'p ^ (? , and is therefore entirely nonrelativistic. A differential formulation that is 
exactly equivalent to the Landau collision integral was given by Franz. ^ 

We turn now to the differential formulation of the relativistic collision integral due to 
Beliaev and Budker.^'^'^ They obtained the expression 



in which 7 = y'l + jup, 7' = + |u'P, and r = 77' — u.u'. (We set c = 1 in this part of 
the paper.) One finds 

£7-'^Ku')=(;J^^(r-u-u'). (4b) 

Notice that r is the relativistic correction factor for the relative velocity between the two 
particles (i.e., for the velocity of one particle in the rest frame of the other). Conversely, 
this relative velocity is given by r~^(r^ — 1)^/^. 

In developing a differential formulation for the collision term based on the Beliaev and 
Budker kernel, it is helpful to work in terms of relativistically covariant quantities. The 
expression 77' 1/ is equal to the space part of a four- tensor W that depends on the four- 
vectors u= (7, u) and u' = (7', u'), 

^2 

W'^{u, u') = _ ^^3^^ ((r^ - l)g'^ - u'u^ - u''u'^ + r{u'u'^ + u''u^)), (5a) 

where g'^^ is the metric tensor, with signature — h++. (r = —Uiu'^ is clearly a four-scalar.) 
The tensor W is symmetric (W-' = W^^), symmetric in u and u', satisfies UiW^ = 0, and 
satisfies W- = 2r^(r^ — 1)~^/^. Likewise ^^{d/du') ■ U is the space part of the four- vector 
V, where 

^'("'"')= (,2_l)3/2 (^^'-"")- (5b) 

If the relativistic differential formulation is to parallel most closely the nonrelativistic 
formulation, then one should find a representation of the form W^^ = W^ip and = —2g^ip, 
where ijj and ip are four-scalars depending on u and u', and and are covariant 
differential operators acting on the variable u. In the nonrelativistic limit, if'j should reduce 
to |v — v'l and if should reduce to |v — v'|~^. It should be possible to transform tjj and 
if to delta functions by a sequence of second-order differential operators. The potentials 
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would be defined as /i = -(I/Stt) /(t/^Z/V) rf^u' and g = -{I/Att) f{tpf/'y')<fiu'; these 
expressions define four-scalars (cf. Rcf. 5). The differential equations satisfied by h and g 
follow immediately from those satisfied by ip and (p. Finally, D would be obtained as the 
space part of —j~^W^h and F as the space part of —'y~^Q'^g. In fact, it will turn out that 
the rclativistic formulation has to be somewhat more complicated, but not fundamentally 
different from the outline just sketched. 

A function of the four-vectors u and u' that is a four-scalar must be a function of 
r = —u.u' alone. The form of the differential operators TC^ and is restricted because 
these should be interior operators on the surface u-^ = — 1 in four-space. In addition, it is 
required that = W'^ and UiW^ = 0. Under those restrictions it is found that the most 
general form of and Q'^ , up to a multiplicative constant, is W-'x = 'C'-'x + Q^Cfl''"' +u^u^)x 
and G^X = + Pu^X- Here, a and (3 are arbitrary constants, and 

r^^X = ig'' + u^u'^W^ + u^u^)^^, + {g'^ + uV)u- 4^ , (6a) 



du 

The spatial part of O^x is i-X and that of lOx is Kx where 



/C'x=(5^'+w^«')|4. (6b) 



dvdv dv dv ' 
Kx = 7-^|, (7b) 



in which v = u/7, and d/dv = 7(/ -|- uu) • d/du. If x is a function of r alone then 
and /C*x = {dx/ dr){ru^ — u"). One is thereby led to the representations 



W'^ = [C^ + g'' + u'u^] v/r2 - 1 

- [C^ - g'^ - u'u^] (r cosh" V - ^/r^ - l) , (8a) 
V' = -2V{r{r^ - 1)^1/^ - cosh" V). (8b) 

These representations for W and V are only suitable for constructing a differential 
formulation of the collision term if the functions that occur on the right-hand sides can be 

reduced to delta functions by some sequence of differential operators. For that purpose the 
contraction L = £■ is needed; in terms of the three-space variables it is 

Lx=(/ + uu):^ + 3u-|^. (9) 

If X is a function of r alone, then Lx = {r^ — l){d'^x/dr^) + ^'''{dx/dr) away from r = 1; at 
r = 1 (or u = u') there may be a singularity. Specifically, it is found that 

L(r(r2 - l)-i/2) = -47r7(5(u - u'), 
[L + l]((r2 - l)-i/2) = -4^75(u - u'), 
L(cosh-V) = 2r(r^ - l)"^/^ 

[L - 3] (r cosh" V - ^/r"^ - l) = 4\/r2 - 1. 
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The explicit form of the differential representation of Eqs. (3) based on the Beliaev and 
Budker collision kerneP follows: The potentials are 

ho = -(l/47r) J (r^ - l)-i/2/(u')/y d^n' , (10a) 

/ii = -(I/Stt) J Vr2-l/(u')/7'rfV, (10b) 

h2 = -(l/327r) J (rcosh^V - \/r2 - l) f{u')/j' d^u' , (10c) 

go = -(1/477) J r{r^ - l)-'/^f{u')/Y d^u', (lOd) 

= -(l/87r) j cosh-V/(u')/7'd^u'. (10c) 

These potentials satisfy the differential equations 

[L + l]ho = f, 

[L-3]h, = ho, (11) 
[L - 3Jft2 = hi, 

Finally one obtains D and F as 

D(u) = -7-i[L + / + uu]/ii+47~^[L-/-uu]/i2, (12a) 
F(u) = -j-'K{go - 2gi). (12b) 

Equations (11-12) together with the definitions, Eqs. (7) and (9), provide the differential 
formulation in the relativistic case. 

In order to proceed further analytically, it is useful to decompose the distribution 
function and the potentials in spherical harmonics, e.g., 

oo n 

f{u,e,^) = Y, E /nm(w)Pr(cose)exp(im<^). (13) 

n=0 m=—n 

Here u = |u| (different from the convention used earlier), 6 is the polar angle, and (j) is 
the azimuthal angle. The equation [L — oi\g = f is equivalent to the system of separated 
equations [L„ - a]gnm = fnm, where 

[L„ - a]y = (1 + n^)g + {2u-^ + 3.)| - (^^^^ + a)y. (14) 

After the change of variable x = sinh^'^u and the change of unknown z = (sinhx)~"j/, then 
the equation [L„ — a.]y = w transforms to [P„ — a^jz = (sinha;)""^, where = a + 1 and 

[Vn - a'']z = ^ + 2(n + l)(cotha;)^ + ((n + 1)^ - a'')z. (15) 

The solution to the homogeneous equation [I>„ — a^]2: = is required in order to construct a 
Green's function for the problem. To obtain this solution we note the following recurrence: 
If Zn-i,a solves [Dn-1 — o^]z = 0, then Zn,a = (sinha;)~-^((i/dx)2;„_i,o solves [P„ — a?]z 
i I. Furthermore, for n = — 1 the homogeneous equation is trivial to solve. However, the 
recurrence breaks down in the case that a is an integer. If a = n, then z„_i,a = 1 solves 



[T>n-i — a'^]z = 0, and differentiation produces the null solution to [Vn — a?]z = 0. The 
recurrence must then be restarted from the general solution to [T>n — 'n?\z = 0, which is 



The integral that occurs here can be expressed in closed form. 

The Green's function allows us to reduce the separated ordinary differential equations 
to quadrature. An important special application for these results is in the treatment of 
collisions off an equilibrium background distribution. Assuming that fb is a stationary 
Maxwellian with density nf, and temperature Ti, and that the energy of the colliding particles 
greatly exceeds Tf,, we obtain 



and Fu = —{mav/Th)Duu- (The other components of D and F vanish.) Here we have put 
the expressions for D and F into dimensional form as in Eqs. (2), Kn is the nth-order Bessel 
function of the second kind, the argument for the Bessel functions is mftC^/Tf,, u^j, = Tf,/mb, 
and Tab = nbqaOb ^og Aab / (^iTTeQm^) . The errors are exponentially small in u/utb- 

To conclude, we have presented a differential formulation for the Beliaev and Budker^ 
relativistic collision integral. This permits the rapid numerical cvahiation of the collision 
term. A decomposition into spherical harmonics is useful in carrying out analytical work. It 
also provides a convenient method for calculating the boundary conditions for the potentials. 

We are grateful to N. J. Fisch for several enlightening discussions. This work was 
supported by DoE contract DE- AC02 76 CHO 3073. 
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